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' By Arthur W. Goldstein 


SUMMARY 

General equations are developed for isentropic, frictionless, 
axisymmetric flow in rotating impellers with blade thickness 
taken into account and with blade forces eliminated in favor of 
the blade-surface function. It is shown that the total energy of 
the gas relative to the rotating coordinate system is dependent 
on the stream function only, and that if the flow upstream of 
the impeller is vortex-free, a velocity potential exists which is a 
function of only the radial and axial distances in the impeller. 

The characteristic equations for supersonic flow are devel- 
oped and used to investigate flows in several configurations in 
order to ascertain the eflect of variations of the boundary condi- 
tions on the internal flow and the work input. Conditions 
varied are prerotation of the gas, blade turning rate, gas velocity 
at the blade tips, blade thickness, and sweep of the leading edge. 

INTRODUCTION 

No method of practical importance has been described in 
the literature for computing the flow in the three-dimen- 
sional region between the blades of a supersonic compressor. 
This deficiency in design technique has been partly met by 
using two-dimensional theory to compute the variations of 
the relative flow from blade to blade, by assuming no radial 
flow components, and by neglecting spanwise variations. 
This procedure is probably satisfactory for designs with a 
ratio of blade-root radius to blade-tip radius close to unity 
and for cylindrical hub and casing shapes. When this ratio 
is small, the spanwise variations are significant. A method 
devised at the NACA Lewis laboratory during 1950 is de- 
scribed herein for computing the magnitude of these varia- 
tions, when blade-to-blade variations are neglected, by 
assuming a large number of blades in the impeller and deal- 
ing with the circumferential average flow. This technique, 
when combined with estimates of blade-to-blade variations, 
is useful for designs with closely spaced blades. 

The usefulness of this approach may possibly be extended 
to supersonic flow in a manner analogous to the method 
described in reference 1 for subsonic flow; the axisymmetric 
solution is first found and then used as a base for computing 
the blade-to-blade variations. The particular method used 
in the blade-to-blade calculation of reference 1 is series 
expansion, which is not applicable to supersonic flows be- 
cause of the existence of solutions with discontinuous 
derivatives. 

In reference 2, the equations for internal flow with zero 


blade thickness are translated into the characteristic equa- 
tions for the derivatives of the stream function, and it is 
shown there that the equations are hyperbolic if the relative 
velocity is supersonic, provided the ratio of the tangential to 
the axial velocity is a given field function. (This condition 
includes the special class of blades containing radii from the 
axis of rotation.) If the tangential velocity is prescribed, 
the meridional velocity component is the decisive factor in 
determining whether the equations are hyperbolic. 

The present report treats the velocity components directly 
after the manner of reference 3 for arbitrary preassigned 
blade shapes. The design problem is solved by finding the 
hub shape that will give a prescribed velocity distribution 
along the casing (or any other selected stream surface). In 
the course of the solution, the tangentially averaged flow is 
completely determined in the interior of the impeller. 

The purpose of this paper, in addition to describing the 
method in a mathematical outline, is to examine the per- 
formance of the technique in computed examples and to find 
whether any peculiarities in the flow arise from the boundary 
conditions imposed, and to develop a background of experi- 
ence for the design of impellers of the type considered by 
examination of the effect of systematic variations of impeller 
design parameters on the internal flow. The parameters 
varied are inlet prerotation, blade turning rate, gas velocity 
distribution at the casing, blade thickness, and sweep of the 
blade leading edge. 

Approximately 160 hours of computing are required to 
work out about 100 points for each of the impellers. This 
time can be reduced by application of the punched-card 
computing technique to part of the procedure. 

DEVELOPMENT OF EQUATIONS FOR AXISYMMETRIC FLOW 

Equations describing the fluid flow are developed in con- 
venient form by use of the equations of motion and con- 
tinuity, the conditions of axial symmetry, and. isentropic 
state changes. The uniformly distributed blade forces, 
which are inconvenient in calculations, are eliminated from 
the equations by substituting the blade surface function, so 
that the condition of integrability is automatically satisfied. 

EQUATION OF CONTINUITY 

The continuity equation for steady flow is 
V(pW)=0 


1 Supersedes NACA TN 2388, “Axisymmetric Supei sonic Flow in Rotating Impellers” by Arthur W. Goldstein, 1951. 
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(Symbols are defined in the appendix.) If axial symmetry 
is assumed, the equivalent in cylindrical coordinates is 

^(pw)+^(pw)=0 

The presence of blades serves to modify these equations if 
the space filled by the blades is taken into account. If the 
impeller contains N blades with thickness t measured in the 
direction of rotation, the ratio of volume available for fluid 
flow to the actual volume of any space element approaches 
(in the limit as the number of blades increases indefinitely 
with Nt fixed) 

A Nt 

r 2irr 

The continuity equation is therefore obtained bj r reducing 
the specific-mass-flow terms pv and pw by the ratio A/r to 
account for the area reduction : 

(p^)+^ (pwA)=o (l) 


Because none of the scalar quantities is a function of 6, 
equation (1) implies the existence of a stream function 
ip(r, 2 ) such that 

\p r =pwA 


\j/ 2 ——pvA >- 

*A«=o J 


( 2 ) 


be shown for frictionless flow) for any particle; energy is 
interchanged among the forms of enthalpy, kinetic energy, 
and potential energy ( — co 2 r 2 / 2) resulting from the centrifugal 
force field. In reference 2, E is shown to be constant for 
any particle for frictionless flow and no heat transfer with a 
finite number of blades in the impeller. 

In order to show that E is constant for any particle in 
axially symmetric flow, the direction of the field force is first 
identified. If the blades are considered to increase indefi- 
nitely in number and to retain their camber-line shape, the 
forces exerted by the blades are merely the pressure forces 
exerted normal to the blade surfaces, which in the limit 
approach uniformly distributed forces in the body of the 
fluid. The relative flow lines are always contained in one 
of the family of blade surfaces, so that the field or blade 
forces must be normal to the velocity. If the scalar product 
of equation (3) is taken with the relative velocity, the result 
is then 

WNE= 0 


which indicates that E does not change along any flow line. 

Elimination of blade force. — The blade force can be 
eliminated from the equation of motion if the equation for 
the blade surfaces is assumed known. The force may then 
be evaluated after the velocity components are known. 

Under the condition of axial symmetry, the equation for 
the family of blade surfaces 

a{r, d, s) = constant 


(subscripts indicate partial derivatives) or 

V<A = P Aj X V= P Aj X W (2a) 

EQUATION OF MOTION 

Substitution of energy for pressure. — The equation of 
motion for steady, nonviscous, isen tropic flow in terms of the 
relative velocity is 

V^—WXVXW+- vp=u 2 R+2WXtt+F 

2 p 

where u?R is the centrifugal force and 2 It 7 XU, the Coriolis 
force. For isentropic conditions, the substitution 

- Vp—Vh 
p 


may be simplified. If the constant is changed, there is a 
shift from one blade to another. Because of axial symmetry, 
this shift is equivalent to rotating the original blade by a 
fixed angle. Therefore, the equation for the surfaces can be 
put into the form 

a = 6+f(r,z) (5) 

Because the relative velocity flow lines are contained in the 
blade surfaces, W is normal to Va as well as to Vi/e The 
relation among the three vectors can therefore be shown 
to be 

V^XVa (6) 

pA 

From the condition of axial symmetry, 


may be made, so that the equation is then 

v(h+^-*~)=WX(vX w +2Q.)+F=WXvXVAF 

( 3 ) 

The quantity 

W 2 V 2 

E^h -\ — ^ 2~ = ^^ — 2 — ' “Ku A ">’) (4) 


is a measure of the total energy of the gas particles in the 
relative system of coordinates and remains constant (as will 


bE 

de 


=rj-VE=0 


Because both VE and Vf are normal to j and W, they 
must be parallel, and consequently E is known to be a 
function of \f/ only (reference 4) ; that is, 

E=E(t) \ 
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Insertion of equations (4), (6), and (7) into equation (3) 
yields (without resolution of the term VX V) 

% VXF ) V “' VXF ) V *+ F 


The blade pressure forces F are normal to the blades and 
therefore parallel to Vet; this relation thus includes only two 
nonparallel vectors Va and V^, so that the components must 
separately satisfy the relations 


F=~(jZ VfVXF) Va=-(rj- W XVX V)Vc. 


and 


dE r _ _ wT7 . 

-rr= t Vet • VX V 

d\p P A 


( 8 ) 


(9) 


Special conditions. — Two special conditions of technical 
importance are considered: In high-speed machines, blades 
with radial elements are desirable in order to eliminate bend- 
ing loads resulting from centrifugal forces. Another fre- 
quently occurring circumstance is that of zero vorticity 
upstream of the impeller. 

For blades with radial elements, the vector i lies in the blade 
surface and therefore 

i .vof=y r =o 

Therefore/ is a function of z alone. From equation (6), 


Then the value of the total energy E relative to the impeller is 


W 2 „2 r 2 T72 

E=h+dE--^-=h+-Y-o 3 (ru+^) 


which is constant upstream of the impeller. Then by equa- 
tion (7), I? remains constant in the impeller, so that the equa- 
tion of motion (9) becomes 

VccVXF=0 (11) 


If the system of coordinates used is r, z, and a for expressing 
equation (11), there results (reference 4) 


A ( r F-VaX*)] r B =[^r (rF-iXVco ] 2 a 


The subscripts indicate which variables are held constant. 
This equation and the fact that the quantities in parentheses 
are functions of r and z only establish the existence of a poten- 
tial function of r and z only (4>(r,z)) such that 

[sl= rF ' ixv “ 

el- 


0 = W ■Va=™-{-vf I -\-wf z (10) 

so that for radial blade elements 

1 f 

u = —wr-jp. (10a) 


The last relation indicates that with r and z constant, 4> is 
constant; this is true regardless of the third independent 
variable, which may be taken as 0. The three equations for 
<f> are then put into vector form 

V<f>=i(rV-VaXk) + Jc(rV-iX'Va) (12) 


This relation provides considerable simplification of the equa- 
tions developed later for solving the flow problem. 

For vortex-free flow upstream of the impeller, equation (3) 
is reduced for conditions there by setting fi=F=VXF=0 
and substituting V for W to get 

v(*+?)=0 

or 

V 2 

h- 1 — — = constant 


The absolute velocity may then be shown to be 

F=V4>-|-r('a-|-cor)Va (13) 

and from 

V • Va= W -Va-f- corj ■ Va= co 

equation (13) gives for the moment of momentum p, 

. . . co — Vb-Va 

p=r(u + a>r) = 

so that 

A. Y 7 ™ . Y 7 dV\ 

(13a) 


(14) 


Tr . (co — Va-V4>) 


(Vet) 2 


Also, zero vorticity gives, for axial symmetry, 

0=VXF— i £(« + «•)+; (!f-£)+t / («■+-*) 

From the first and last terms and axial symmetry, 
ur + cor 2 = constant 


The potential can be used as an alternative to the stream 
function in solving for the flow. Equation (13) indicates 
that the projection of V4> on the blade surface is equal to 
the projection of the absolute velocity on the blade surface. 
The equation for ^ and that for 4> are obtained in the following 
manner: The velocity V is eliminated from equation (9) by 
means of equation (6) and there is obtained 
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T “ • vx Gx W-XV„)+2„/,=-fA (g 

which is also applicable in the general case with nonzero 
vorticity upstream of the impeller. 

In order to obtain an equation for <h, equation (13a) is 
substituted into equation (2a) 


/ w 2 , ruwi 2 \ , / ruv f z wv\ . / ruw f T wv\ . 

( 1 “2 4 - — 3 '- ) v I + r A(,w 2 J l t+2wvi, l + »*]„)+ 

\ CL CL / CL 

K i+ ^"S) + z[^ ( ' 4/r)+w, ^ (A/7,) ] =o (i7) 


Then 


Vp = pArve XF= pArVO XV’t>+ pAprVd XV a 
0=VXVp=VX[pArVOX(V3>+p.Vct)] 


and for the equation of motion 


( 1+ r 2 /z 2 ) W T + r 2 Jzfr( V r — w z ) — ( 1 + f *f 2 ) V z + 


Although obtaining a solution in terms of p is feasible 
whether dE/dip is zero or not, when dE/di/ = 0, <f> may prove 
more convenient. In regions where IT 2 /® 2 is nearly 1.0, 
because |IT| is a double- valued function of the stream- 
function derivative p| IT|, it will be uncertain whether the flow 
is subsonic or supersonic. A solution obtained from a net- 
work would require a very fine spacing in regions where 
|IT| —a. 

GENERAL EQUATIONS IN CYLINDRICAL COORDINATES 

In general, where f T X0 and dE/dpAO, the equations for the 
internal flow are solved simultaneously for the two velocity 
components v and w, when u has been eliminated from the 
equations by using equation (10). The equations are trans- 
lated from their vector form to scalar notation by utilizing 
cylindrical coordinates. Then a further transformation is 
employed with characteristic quantities as the independent 
variables. 

The continuity equation (1) is first written in vector form 
(the condition of axial symmetry being used) and multiplied 
by r/pA 

-^r V -ip - w)= V • W+ IT ■ ~ +4 W ■ V(A/r) = 0 
pA \ r / p A 





)] 


— 2/x 



. dE . 
pA ~dp = 0 


(18) 


USE OF CHARACTERISTIC VARIABLES 


Statement of basic equations.- — The method of reference 3 
is applied to equations (17) and (18). The object is to find 
two alternative equations with new variables £ and ??, instead 
of 2 and r, such that each equation (corresponding to equa- 
tions (17) and (18)) contains derivatives with respect to only 
one of the independent variables. Two additional equations 
are required to find the position coordinates z and r in terms 
of the characteristic coordinates £ and r\. There are two 
characteristic slopes dr/dz, and f_, the first corresponding 
to £ varying and -q constant, the second to £ constant and 17 
varying. If f represents either of these slopes and a repre- 
sents either of the coordinates corresponding to the approp- 
priate f, the results are 

r a =$z„ ( 19 ) 

r =4 r\Vaf± Vr 2 (V a)\M 2 — 1 j] (20) 

(the upper sign of ± refers to f + , the lower, to f_) and 


For isentropic flow, 


V p Vh 


Introduction of E by equation (4) and use of IT- VE=Q 
produce 


V-IT- 


IT- VlT 2 /2^o) 2 7 , «_ L r ^ 


- wZ fr)=° 


a 2 n a 2 ' A\ dr 
This expression finally reduces to 

--?(£), -S(f) > + 

. <15 > 

In cylindrical coordinates, equation (9) becomes 

dE 


Qr Vz~\~j Cpz J zPt — pA 


dip 


(16) 


If equation (14) is now substituted into equations (15) and 
(16), there results for continuity 


where 

B=l+r 2 f 2 -^ 


J WaAC v Lz „— 0 


r 2 (Va) 2 


( 21 ) 


M 


J 


2 __(ti 2 +W 2 + W 2 ) 




G = (M 2 — l)r 2 /J z +^r 2 (Va) 2 ±Vr 2 (Va) 2 (M 2 - 1) 

L^-(?-\fjA+l+r 2 J 2 )D l + 

Vv, , . A w 2 1 ruw , \ ."] „ 

D i =™(w 2 f zz + 2w v f TZ + V 2 Jrr^) + 

v / Ar 2 u 2 \ . vr d(A/r ) , wr d(A/r) 
r \ ‘ ad a 2 )-^ A dr + A dz 

A-r 2 / z 2 [40+ W ^)]- 2( u+ W r)/ 2 -p^^f 


Y (21a) 
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Equation (20) for f shows that the system is hyperbolic if 

M 2 > 1 . 

In the homogeneous equations (D l =D 2 —0), the coefficients 
L + and (corresponding to f + and f_) are zero; conse- 
quently, in the hodograph plane the slope of the characteris- 
tics would be given by — plane flow, JJG is a 

function of v and w only, and this equation may be integrated 
and the hodograph solutions tabulated in advance of applica- 
tion to any particular problem that may arise. In the pres- 
ent circumstances, however, J/6 also involves the position 
coordinates, and 0, Z? 2 A correspondingly easy 
method of solution by simple waves is therefore unavailable. 



Figure t. — Characteristic and cylindrical coordinates. 


Instead, the four (<7 = £ or 17) equations (19) and (21) arc inte- 
grated simultaneously in a step-by-step procedure. 

Method of computation. — In figure 1, conditions arc sup- 
posedly known at P, and P 2 and conditions are to be computed 
at P 3 . Then equation (19) in difference form is 

r 3 —r l =(z i —z 0f + (22) 

n~r 2 ={z 3 —zd)^ (22a) 

This system may be solved for r 3 and z 3 . Then similarly for 
equation (21), the difference form is 

J f V) 3 —w L ) + G + (v 3 —Vi) +£+( s 3 — s,) =0 (23) 

J(' w 3 — w 2 ) + 6?_ (» 3 — v 2 ) + X_ ( s 3 — z 2 ) = 0 (23 a) 

The subscript + indicates that 7? is constant and the subscript 
_ indicates that £ is constant. In equations (22) and (23) 
the coefficients f, J, G, L are averages between values for the 
points Pi and P 3 . For the first approximation, only the values 
at Pi are known and these values are used. 

When approximate velocity components and position coor- 
dinates are thus determined for P 3 , the coefficients f , J, G, L 


may be computed there and averaged for a better approxima- 
tion of the position coordinates and the velocity components 
at P 3 . This process is repeated as often as necessary to 
arrive at an answer of satisfactory accuracy. The smaller 
the interval, the quicker the convergence to the solution. 
In each stage of the process, if dEjdij/^0, then ^ 3 must be 
computed from 

^3 — ('Prt + 'Pz)dz = (\p T + &/ f)dr - 


From the value of \J/ 3 , dE/dip at P 3 may be found.. ; 



Figure 2. — Relation between shape of boundary and boundary conditions. 

BOUNDARY CONDITIONS 
GENERAL 


The discussion in the previous section implies that boundary 
conditions are known on a curve C containing points such as 
Pi and P 2 , and that the solution is then extended into a region 
con taming points such as P 3 . The shape of the blades and 
E 0/0 are assumed known in that region. On the curve C, 
there are known in addition the values for w and v, from 
which all other variables (such as dE/d\f/, \j/, u, p, a 2 , and so 
forth) may be computed. Suppose the initial curve is such 
that the characteristic curves intersect it at two points, such 
as A and D or A and E in figure 2. Then by equation (21), 
there is a relation between the variables at A and D that does 
not permit the independent specification of the variables 
v and w at A and D. Similarly if two characteristics passing 
through two boundary points reach a third boundary point, 
then another degree of freedom is lost. In figure 2, if two' 
variables are prescribed on A B and EF, then one variable may 
be prescribed on BC and DE, and none on CD. 

TIP CONDITIONS 

The supersonic impeller has a high pressure ratio per stage 
and a large air-handling capacity per unit of frontal area. 
In order to maintain the large capacity, the tip radius at 
the entrance should be at least as large as at the exit. In 
order to provide a maximum work output with a minimum 
turning, the exit tip radius should be at least as large as the 
entrance tip radius. A reasonable shape for the tip stream- 
line is therefore a surface of constant radius, although smaller 




6 


REPORT 1083— NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


radii between the entrance and the exit might be useful in 
order to affect the internal flow. If the casing shape is se- 
lected, then, as it is also a streamline, the ratio vjw is 
specified. Selection of the velocity or of a velocity component 
remains. This selection may be made by consideration of 
the blade-to-blade velocity variations and by using some 
approximate method such as that outlined in the following 
paragraphs. 

In three-dimensional motion (finite number of blades), 
Kelvin’s theorem states that if the absolute velocity is irro- 
tational upstream of the impeller, it will remain so inside the 
impeller ; that is, 


from the axisymmetric solution, and further that (W'-g),= 
\W/\, then 

|ir,1-|ir/|-g + 0+2„) T A f ]?^ (24a) 

if the relation 

v /\W\=i-q 

is utilized. In this equation all the velocity components 
without subscripts are obtained from the axisymmetriG 
solution. If 

\W' S \ + \W' P \ = 2\W\ 


VXV'=VXW' + 2cok = 0 

where the prime indicates velocities of the nonsymmetric 
flow. Then 


may be reduced to 


V^-VX K' = 0 


A (W'-g)=r (u' + 2a,r)g_-i 


(24) 


where bu'/ds is the arc-length partial derivative of u' 
along the direction of #. This equation is now integrated 
from the pressure surface of one blade to the suction surface 
of the next in the direction r= constant, 2 = constant. The 
variable 6 increases by the amount 2-wAjNr, so that 

(W'>a) s -(W'-q) P ==rj e de+cl-ij^ u'd,6- 2c ,rg_-i 

If it is now assumed that the values of du'fd-s and of u' 
averaged across the channel are equal to the values obtained 


A B C 



(a) 



(a) Conditions on B'B and BC unrelated, 
(b) Conditions on B'B related to those on BC. 
Figure 3. — Boundary conditions for impeller. 


or some other connection between |K T/ S |, |TT'p|, |W|, and 
the impeller geometry is assumed, the blade surface veloci- 
ties at the casing may be obtained with that relation and 
equation (24a). The distribution of |TI 7 | at the casing may 
therefore be selected on the basis of obtaining a desired 
blade surface velocity distribution or distribution of u. 

Suppose that, in some manner, all the variables have 
been selected at the casing (curve BC in fig. 3(a)). The 
flow is then uniquely determined in the region BPjC enclosed 
by the casing and two characteristic lines. If the flow is 
also prescribed at the leading-edge surface BB', then the 
flow is determined in the extended region B'BCP 1 P 2 . When 
the characteristic line BP 3 falls inside this region (fig. 3(a)), 
no difficulty results. If, however, the characteristic B'P 4 
intersects the casing inside the blade region (fig. 3(b)), then 
the velocity may not be prescribed on the curve BP 4 (refer 
to discussion of fig. 2), or alternatively, only one of the 
velocity components may be prescribed on BB' (BB' is not 
a streamline) . Similar considerations are valid on the curve 
CC' at the exit. 

The selection of values for the variables v and w on the 
entrance edge BB' must be made by joining the solution 
for the flow in the region ABB' A' with that for the flow in 
the blade region BCC'B'. Because the blades may deflect 
the air at BB' discontinuously and because the effect of 
thickness may be felt by the air in a discontinuous manner, 
the discontinuities in the flow should be evaluated. A 
method for making this estimate follows: 

The angle is defined as the angle between the meridional 
projection of the blade leading edge and the z-axis (fig. 
3(a)). The conditions outside the blade region are desig- 
nated by subscript zero, those inside the blade region, by 
no subscript. Then across the surface BB' the equation 
of continuity yields 

p 0 (w 0 sin w 0 cos p)=(l— ) p( w /? — v cos /3) (25) 

The discontinuity in blade thickness can cause no force 
in the entrance surface BB', and the blades can exert no 
force in the blade surface. Therefore no force and no 
discontinuity in momentum exist in the direction of the 
curves representing the intersection of these two surfaces; 
that is, 
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(W cos fi — k sin (3)XVa= 0 
or 

r(u—Uo)(f r sin /3 +/* cos /3)=(» — w 0 )sin p + (w— w 0 )cos /3 (26) 

(27) 


(28) 


There are then available equations (25), (26), (27), and 
(28) for finding p, u, v, and w. The initial values of f r and J z 
may be chosen from these equations to obtam desired rela- 
tions between the velocity components outside and inside the 
blade regions. 

APPLICATIONS OF METHOD OF CHARACTERISTICS 

In the examples computed in this section, the conditions 
for which solutions are obtained vary. The purpose of these 
solutions is to find the effect of some variations of prescribed 
conditions on the resultant flow. For all the examples, 
dE/d\f/=0; the flow is vortex-free upstream of the impeller; 
at the leading and trailing edges the radial component v of the 
velocity is zero; and the blades contain radial elements 
(J r — 0). The relative Mach number at the conjunction of 
the leading edge of the blade and the impeller casing is 1.94. 
Except in example VI, the value of the blade angle with 
respect to the axial direction is zero at the exit (j t — 0). The 
other conditions are summarized in the following table: 


There is also available 

u=—r(wf z + vf T ) 

If the process is assumed isentropic, then 


_P 

Po 


Kr-o- 


iiy- if 2 \ 


7-1 


2Ao / 


is obtained by utilizing the fact that 


w , , UV , , W 2 <j?r 2 

E=h 0 +— __ 


and that 


r 0 —r 


Example 

Entrance 



Casing 

velocity 

Turning 

rate 

Sweep- 

back 

(deg) 

Blade 

thick- 

ness 

Prerota- 

tion, 

(w+wr)r 

Axial 

inlet 

velocity 

Tip 

speed 

Blade 
angle 
at tip 
(deg) 

I 

0 

0.85 ai 

1.68 aj 

-61.8 

(*) 

0 to 

maximum 
to 0 

0 

0 

II 

—0. 3 airt 

0. 835 ai 

1. 44 ai 

-64.4 

1.94 ai 

0 to 

maximum 

toO 

45 

' 0 

III 

—0.3 air* 

0. 835 ai 

1. 44 ai 

-64.4 

1.94 ai 

Twice 

example 

II 

45 

0 

IV 

—0.3 am 

0.835 ai 

1. 44 ai 

-64.4 

(»> 

Twice 

example 

II 

45 

0 

V 

—0. 3 din 

0. 835 ai 

1. 44 ai 

-64.4 

oo 

Twice 

example 

II 

45 

w 

VI 

—0.3 am 

0. 835 tti 

1.44 ai 

-64.4 

(») 

0 to 

maximum 

10 

0 


* Determined by conditions at leading edge to z/rt= 0.356; held constant at 2.1 a from 
z/r<= 0.356 to 1.0. 

b Constant at 1.94 a to z/n= 0.25; decreased to 1.4 a (=1.6 aO at exit. 

o Zero from z=0 to z/rt= 0.25; from z/r t = 0.25 to z/r t = 0.5 the correction for blade thickness is 
A=r exp |~~~ I^I—cosStt ^^—0.25^ JJ 
d Determined by leading edge; initial value, 1.94 a. 


Example I.- — The results of the computation for the first 
example are summarized in figure 4 (a), which shows stream- 
lines, contours of constant Mach number and constant energy 
level, and some of the characteristics of the ^-family, which 
intersect the leading edge. Because of the intersection of 
the characteristics from the leading edge with the tips of 
the blades, there is no control over the Mach number, as the 
blade shape was prescribed by 


d6 dj , nA r /i I 

r -j — = — r ~r~= — 1.04 — 1 1+cos 
dz dz r t \ 



The Mach number at the casing rises to a maximum of 2.1 at 
s/r,= 0.356, and is then assigned a constant value. By equa- 
tion (4), the stagnation enthalpy H ( = h J r~V 2 l 2), which 
measures the energy of the gas, increases with z at the rate 

bH du 

bz bz 


At the casing »=0, and therefore, from equation (10a) and 
W 2 —u 2 -\-w 2 , 


did t or 

bz \W\ 



rw 3 d 2 J\ 
\W\dz 2 ) 


This expression shows that the rate of work input is increased 
by increasing the rate of blade turning and decreased by 
increasing W. Figure 4 indicates that the rate of turning of 
the blade is Loo small in the beginning of the impeller and 
that as a consequence nearly all the work input is concen- 
trated in the back portion of the impeller. This condition 
is aggravated by the acceleration of the gas in the initial 
portion (w]IF| 2 <(0 because u<C0) . Increasing the turning 
rate, however, will open the flow area and cause \W\ to 
increase unless the area is controlled by blade thickness. 

A great disparity of work exists between the root and the 
tip; the values for AH/ai 2 are 0.65 and 2.5 at the exit. For 
isentropic compression, the corresponding ratios of outlet to 
inlet stagnation pressures are 2.48 and 9.0 for the root and 
the tip, respectively. The Mach number at the root at the 
exit ( z=r , , r=0.507 r ,) is in the neighborhood of 2.25, and 
a maximum value of about 2.5 is reached at z=r t , r= 0.70 r t . 
Subsonic Mach numbers are reached in the region 
0O/r,<0.3. The solution could not be extended into this 
region. 

If the velocity were caused to decrease at the exit ( z=r , ) 
by decreasing the velocity at the casing ( r=r t ) in the rear 
part of the impeller (0.356<2/r ( <1.0), the mass-flow density 
at the exit would increase with the decrease in velocity and 
the required flow area would decrease. Termination of the 
hub would be expected, not at r=0.507 r t , but at some larger 
value where the work output would be higher. Another 
method of reducing the blade height at the exit and thereby 
reducing the work variation there would consist in relaxing 


219512—52 2 
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AH/a* z 



(e) Example III. 

Figure 4. — Continued. Flow characteristics. 


the condition of v=0 and utilizing some other relations that 
would correspond to lower Mach numbers. Two variables 
or relations could also be prescribed at the exit, but conditions 
on the blade tip (r=r t ) could not be specified for some region 
near the back of the impeller. In this respect the entrance 
and exit conditions would be similar. 

Example II: Prerotation of gas. — In order to maintain 
supersonic flows throughout the impellers so that the flow 
could be computed, the relative entrance Mach number at 
the blade root for example II was increased by assuming an 
initial counterrotation. The rotative speed had to be reduced 
from 1.68 ay to 1.44 ay to maintain the work output, the 
relative Mach number at the leading-edge tip, approximately 
the same blade angle there, and the same axial velocity com- 
ponent at the inlet. This reduction had a further beneficial 
result of reducing the absolute Mach number at the exit for 
a given relative velocity axially directed. The condition of 
2 = constant for the leading edge cannot be maintained for 
shockless entry with radial blade elements and a counter- 
rotating free vortex entering the blade system. The sweep 
angle /3 chosen for this example was 45°, enough to permit 
free selection of TT= 1.94 a over the entire blade tip, because 


the characteristic curves from the leading edge do not inter- 
fere with boundary conditions on the casing. The reduction 
in velocity at the casing reduces the blade height at the exit 
and the work-output variation there. 

In order to examine deviations from an example more 
nearly representative of a good design, the velocity at the 
blade tip for example II was not maintained the same as in 
example I, and consequently the differences between examples 
I and II are not ascribable solely to inlet prerotation. Com- 
parison of the effect on rotative speed (for the same turning 
at the tip and the same work output) and on diffusion at the 
blade roots is, however, valid. 

The results of the computations are shown in figure 4 (b). 
In this example, the turning rate is predetermined in the 
region — 0.5 < O/r,<(0 by the velocity distribution upstream 
of the leading edge, the condition of radial blade elements, 
and the slope of the leading edge. The turning rate is con- 
tinuous at z— 0, and satisfies the equation 

i (^) =3 - 008 ~ 2 - 298(3s-1)2 ( 0 <«/ r ,< l / 3 ) 

= 3.008 [1-1/4 (3 2 -l) 2 ] (l/3<s/r i <l) 
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Compared with the first example, the work-input rate at the 
blade tip is greater in the initial section; however, the work 
distribution is still unsatisfactory. The exit blade height is 
also much less, with the result that the work input is much 
more nearly equalized; at the root the value of A H/a? is 1.8, 
whereas at the tip the value is 2.5. These values correspond 
to stagnation pressure ratios of 5.2 at the root and 8.35 at the 
tip. This variation could be reduced by sweeping back the 
trailing edge to provide space for additional work input at 
the hub streamline. This sweepback is more feasible with 
the smaller blade height and work variation of the second 
example than for the first example. 

Example HI : Effect of turning rate. — In order to reduce the 
size of the impeller by increasing the turning rate, example 
III was computed with twice the blade turning rate of exam- 
ple II for s>0 and with the same over-all work output. One 
result is the greater diffusion on the hub, where a minimum 
Mach number of 1.2 is reached instead of 1.4 as in the pre- 
vious example (see fig. 4 (c)). The increased diffusion 
probably results from the increase in hub curvature, which 
in turn is caused by the shorter axial distances required for 
the changes. Another result is the increased height of the 


blade at the exit, and an increase in the exit Mach number 
from 2.3 to 2.7. The increased streamline curvature also 
accounts for the higher exit velocity and this higher velocity 
in turn explains the larger blade height. The characteristic 
line through the casing at the exit (r=r t , 2=0.5 r t ) intersects 
the hub at 2=0.155 r t . The characteristics through the re- 
maining portion of the hub (z^> 0.155 r t ) are affected by con- 
ditions at the exit (2=0.5 r t ) where the velocities increase 
with decreasing radius. The result is a large increase in 
velocity at the blade hub in the region 2>0.155 r t . 

Example IV : Effect of diffusion at casing.- — In the next 
example, all conditions remained the same as in the third 
example except that the prescribed velocity at the casing was 
varied smoothly from the value of 1.94 a at 2=0.25 r t to 
1.40 a at 2=0.5 r,. The desired reduction has been obtained 
in blade height and work-output variation, but the discharge 
velocity has a large value of 2.6 a (fig. 4(d)). It may there- 
fore be concluded that the alternative of directly prescribing 
a condition with low relative velocities at the exit would have 
been more effective than prescribing conditions at the blade 
tip for reducing the blade height at the exit. 



(d) Example IV. 

Figure 4. — Continued. Flow characteristics. 
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A new phenomenon has also arisen in that the character- 
istics passing through s=0.5r„ r=0.92r, and 2=0. 5r ( , 
7-=0.893r, intersect at z=0.397r,, r=0.832r,. At this 
point, the Mach number changes discontinuously from 2.13 
to 2.43 and the streamline inclination to the axis changes 
from 14.9° to 7.4°. This single point does not represent the 
full extent of this difficulty. Actually, there is a region 
bounded by a curve in which two solutions are obtained; 
each solution is determined by a different region of the 
boundary. Consequently, the boundary conditions are 
incompatible with obtaining a single-valued continuous 
solution. This situation is analogous to the Prandtl-Meyer 
flow, where the characteristics intersect at the corner. If a 
computation were started for a simple two-dimensional 
supersonic flow from a streamline shape and from the velocity 
distributed thereon, which were exactly the same as on one of 
the streamlines of a Prandtl-Meyer flow around a corner, 
then, when the solution was extended toward the corner, the 
characteristics would intersect at the corner in the manner 
of figure 4 (d) . The intersection on figure 4 (d) can be taken 
to represent a physical boundary to the fluid in that the 
streamline passing through that point would be taken for 
the hub shape. Only 90 percent of the flow (iA=0.9) could 
be accommodated in the impeller of example IV if continuous 
flow with the imposed boundary conditions were required. 
An obvious method of increasing the flow is to increase the 
slope of the characteristic lines at the exit for r<0.9r,. This 


increase could be accomplished by decreasing the relative 
velocities there and discarding the condition v=0 to obtain 
a shorter blade height. 

Example V : Increasing blade thickness. — In none of the 
examples computed so far was blade thickness considered. 
An obvious effect of increasing blade thickness is to decrease 
the velocity in a region where there is no independent con- 
trol of the velocity (such as the region 2<(0.35r, in example I, 
fig. 4 (a)). This phenomenon provides the degree of freedom 
required in such regions to control the relative velocity |W| 
as well as the surface velocities |W,| and \W P \. In regions 
where the velocity \W\ at the blade tips is under control (the 
entire tip sections for examples II, III, IV), altogether dif- 
ferent effects may be expected. The effect of thickness 
under these circumstances was investigated by computing 
another example with a cosine-like bump in the blade thick- 
ness starting at s=0.25r,. (See the previously presented 
table of boundary conditions.) The maximum thickness 
was reached at 2 = 0 . 375 ?',. At r=r t the thickness is 20 
percent of the blade spacing, and at the hub the value is 25 
percent. The results of the computation are shown in figure 
4 (e). The effect of the bump is not felt upstream of the 
characteristic with positive slope passing through 2=0.25?',, 
r=r,. The most obvious effect is a decrease of flow Mach 
number in the region of increasing thickness (and inward 
from the casing, where the velocity is not preassigned), as 
might be expected from one-dimensional theory. At z= 



(e) Example V. 

Figure 4. — Continued. Flow characteristics. 
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0.2 5r t and z=0.375r,, the change in blade height is about lOper- 
cent as compared with example IV. Another effect of 
thickness is to be seen in the waviness of the hub streamline 
in the region 0.15<Cz/^(<^0.25 where the thickness has not 
yet started to increase. This bump may be traced in other 
streamlines along the characteristic running through the 
point of maximum thickness at the casing (r=r t , z= 0.375 r t ). 
The effect of the bump on the hub is therefore to increase the 
local velocity by curvature in order to hold the required 
velocity at the casing; this action from hub to casing takes 
place along the characteristic line as simple waves in two 
dimensions, and the velocity is seen to have increased all 
along this characteristic line. The intersection of the char- 
acteristics and the double-valued solution occur near the 
exit over a considerably larger portion of the flow than in 
the previous example. In this example also, if a lower ve- 
locity were prescribed at the exit, a unique solution would 


be obtained. It would also be feasible to carry out this 
idea and at the same time to prescribe another condition 
(such as »=0) provided that one were willing to lose control 
of the condition at the back portion of the casing ( r—r t ) in 
a manner analogous to that used at the front portion in 
example I. 

Example VI: Effect of leading-edge sweep, — It is of more 
technical interest, however, to carry out this scheme at the 
entrance, for then the leading-edge sweep may be reduced 
from 50 percent of the axial depth of the machine to consider- 
ably less. A computation of the inlet flow conditions for 
such an impeller is shown in figure 4 (f) , where the conditions 
at the leading edge are the same as for example III. The 
blade shape is different, with the second derivative of angle 
d 2 d/dz 2 = — d 2 _f(z)/dz 2 constant. Also, the velocity at r— 1.0 
cannot be prescribed. Only in the region enclosed by the 
casing and leading-edge characteristics was the computation 



(f) Example VI. 

Figure 4. — Concluded. Flow characteristics. 
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completed. The result (fig. 4 (f)) shows a large increase in 
Mach number, especially at the blade tip. Because of this 
increase, the blade turning accomplishes no work hi the 
beginning of the impeller, but shows a region of energy 
abstraction from the air; the initial negative rotational com- 
ponent becomes more negative as the air tons axially and 
the velocity increases as in example I. The situation is 
quite different in the back portion where the blade angles 
approach zero and the tangential component does likewise; 
all the work input is concentrated in the back portion of this 
impeller as in the other examples. Postponement of work 
input to the back portion of the impeller is to be expected in 
any case where the velocity is increasing as the blades turn 
toward the axial direction, or where the velocity is decreasing 
as the blades turn past the axial direction. Internal flow 
and work distribution could be considerably improved in 
this example by increasing blade thickness in the region of 
large velocity increase. 

CONCLUSIONS 

General equations developed for isentropic, axisymmetric 
flows in rotating impellers with any given blade shapes 
showed that the blade-force term can be eliminated and 
solved for, if desired, after the flow solution is completed. 
Further, if the flow upstream of the impeller is vortex-free, 
then a potential function of only radial and axial distances 
exists, which may be used to find the projection of the 
absolute velocity on the blade surface. 

For supersonic flow, the method of characteristics was 
used to compute the flow for several examples with radial 
blade elements and vortex-free inflow. These examples 
show that 

(1) With zero radial inflow, the influence of the leading- 
edge and the blade turning on the casing velocity in regions 
where this velocity cannot be prescribed in advance causes 
an increase in velocity and a reduction of blade loading and 
work input in the initial portion of the compressor. These 
effects can be compensated for to some extent by increasing 
blade thickness. 

(2) Sufficient leading-edge sweepback adds a degree of 
freedom to the boundary conditions at the blade tips and 


therefore permits a more constant or a decreasing velocity 
inside the impeller. 

(3) Reduction in the exit velocity results hr a shorter blade 
height and a more uniform distribution of work input at the 
exit. 

(4) Increase in the rate of turning causes more rapid hub 
curvatures and consequently larger internal velocity gra- 
dients. Diffusion on the hub is therefore greater in the 
initial part of the impeller, and a lower velocity minimum is 
reached there, whereas the gas velocity at the hub near the 
exit reaches larger values. This situation requires larger 
blade heights at the exit with a larger exit velocity at the 
root and greater variation of work output from root to tip. 
Diffusion of the gas velocity at the tip will moderate this 
effect, as will also a relaxation of the condition of zero radial 
velocity at the exit and the substitution of a condition giving 
lower exit velocities. 

(5) Not every set of boundary conditions results in a 
continuous, single-valued solution in the region of interest. 
In the examples considered, a relaxation of the condition of 
zero radial velocity at the exit would have permitted another 
condition giving lower velocities there. This condition 
would have eliminated the region of the double-valued 
solution and given a shorter blade and a more nearly 
equalized work output. 

(6) The effect of variations in blade thickness in regions 
where velocity is also prescribed at the blade tip is to change 
the velocity as in one-dimensional flow but not by the same 
amount, for there is some compensation by variations in 
blade height. The prescribed velocity and thickness require 
hub curvatures that cause velocity changes to be projected 
like simple waves on characteristic curves. It is therefore 
possible to find regions of increased velocity when the blade 
thickness is increased. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 24, 1961. 


APPENDIX 


SYMBOLS 


The following symbols are used in this report: 

A r — Nt./2ir 

a sonic velocity of gas 

a\ sonic velocity of gas at conjunction of leading edge 

and impeller case 

a * critical speed of gas upstream of impeller 
B an abbreviation (see equations (21a)) 

D,,Do abbreviations (see equations (21a)) 


E 


F 

i 


sum of stagnation enthalpy and potential energy 
of gas (specific) in moving coordinates 


(~ A+ v- 


■) 


distributed blade force per unit mass of fluid (a 
vector) 

part of blade surface function dependent on r and z 
only ( a=e+f(r,z )) 
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G an abbreviation (see equations (21a)) 

H absolute stagnation enthalpy (h-\-V 1 2 l2) 

AH work input to gas (rise in absolute stagnation 
enthalpy) 

h static enthalpy of gas (specific) 

1 unit vector parallel to curve 0= constant, z= 

constant (i=Vr) 

J an abbreviation (see equations (21a)) 
j unit vector parallel to curve 2 = constant, r= constant 

(j=rV0) 

k unit vector parallel to curve r=constant, d— 

constant (k—Vz) 

L an abbreviation (see equations (21a)) 

M local relative Mach number 

N number of blades 

p static pressure of gas 

2 unit vector parallel to W 

R vector of r 

r normal distance from axis of rotation 

r, distance from axis of rotation to blade tip at leading 
edge 

ds arc-length element parallel to W 
t blade thickness measured in direction of rotation 

u rotational component of relative gas velocity 
V absolute gas velocity (a vector) [iv -\rj (u + &>r) + kw] 
v radial component of gas velocity 

W relative gas velocity (a vector) (iv+ju-\-kw) 


w axial component of gas velocity 
2 distance measured parallel to axis of rotation 

a function for blade surface (a (r, 6, z)= constant) 

0 angle between meridional projection of leading edge 

and axis of rotation 

7 ratio of specific heats 

f either f + or in general form of equations 

slope of curve for J variable, n constant 
slope of curve for ij variable, £ constant 
v one of characteristic coordinates for gas flow in 

impeller 

8 meridional angle (for cylindrical coordinates) 

yu absolute moment of momentum (r (u+cor)) 

£ one of characteristic coordinates for gas flow in 

impeller 
p gas density 

a either i? or £ in general form of equations 

<i> velocity-potential function of r and s only that exists 

for initially irrotational absolute motion 
tp stream function (Vi/'= pAjX V) 

vector form of w 

u> angular veloci ty of impeller 

Subscripts: 

0 upstream of impeller 

p pressure surface of blade at tip 

$ suction surface of blade at tip 

Subscripts r, d, z, a, £, -q indicate partial differentiation. 
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